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! 1 Abstract 

The geodesies for classical particles in the gravitational field of a Schwarzschild 
black hole immersed in a vacuum of non-zero A are found in their functional 
form. In particular, when r is a function of 4», the functional form is in terms 
of quadrupoly periodic double theta functions. The turnig points are branch 
, places on a one complex dimensional dissected Riemann surface. The four peri- 

■ ods are found by integrating across the the period loops. It is shown that when 

, A is equal to zero, the quadrupoly periodic functions for the geodesies reduce 

to the functional form of doubly periodic elliptic functions. 

o 

O 

"o ! 2 Introduction 

5^ \ Steven Weinberg, [Weinberg, l],refered to the cosmological constant problem as 

bJ) 1 crisis in physics, due to the discepency of some 120 orders of magnitude be- 

tween the quantum field theory calculation and the cosmological calculation of 
the energy density in vacuum[Carroll,Ostlie,2]. Some progress is being made in 
/\ ' astronomy with regards to this problem. The findings from the Type la Super- 

nova data,[Filippenko,3], has created a new cosmology which is as revolutionary 
as Hubble's finding's of an expanding universe. Their observations suggest that 
the expansion of the universe is accellerating 

Presently there several investigations in progress [Efstathiou, 4] which at- 
tempt to measure cosmic acceleration, also known as Dark Energy. If this 
is to be modelled by a cosmic force, F = | A r m c 2 where c is the speed of 
light then the observations are in favor of a positive cosmological constant, A. 
This constant will be a fundamental constant of nature, and the cosmic force 
will be a fifth fundamental force of nature. Since the force is proportional to 
r, the method which probes the greatest distances will work the best, such as 
looking at anisotropies in the cosmic background radiation [Efstathiou,4] . The 
effects of the cosmological constant are neglible on the scale of the solar sys- 
tem[Neupane,5][Cardona,6]. It is vital to find independent experiments to mear- 
sure cosmic accelleration[Boughn, Crittenden, 7]. The equations in this paper 
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may be tested on the scale of our galaxy, thus offering an independent measure- 
ment of the cosmological constant, in relatively well known setting. For example 
the work done with galactic velocity rotation curves [Whitchouse, 8] To do this 
consider the nucleus of the galaxy to be a black hole of radius r g = 2GM/c 2 , 
where M is the estimated mass of the galatic nucleus, and let the orbiting par- 
ticle be a star or globular cluster. With such an experiment, one will have 
to make a reasonable estimate of the dark matter present. The mounting ev- 
idence in support of a cosmologicl constant makes it is necessary to know the 
mathematical relations and the geometry associated with an objects motion in 
space of intrinsic curvature and gravity. The dcvclopcmcnt of the mathematical 
formula that detemine the geodesies where formulated mostly by Riemann and 
Weierstrass, at least the general case[Baker,9,10]. 

To calculate the geodesies in the centrally symmetric gravitational field of a 
Scharzschild black hole in a universe with a cosmological constant, A ^ 0, the 
Hamilton-Jacobi method is applied to the metric, [Witten, 11], [Hawking, 12], [Rindler, 13] 

ds 2 = g tt dt 2 + g rr dr 2 + d<p 2 + g eg d6 2 (1) 
9tt = (l-^Ar 2 - r g /r) c 2 

9rr = -(1 - i(T A r 2 - rg/r)- 1 

g<t><t> = -sin 2 9 r 2 
9ee = -r 2 

Where u = ± 1 depending on whether or not the universe is de Sitter or anti-dc 
Sitter, respectively. A is the cosmological constant and was first introduced by 
Einstein [Filippcnko, 3] and r g was introduced by Schwarzschild. The functional 
form of the geodesies when r g = have been studied, they are given in terms 
of genus zero circular functions. When A = the geodesies are given by genus 
one elliptic functions [Bartlett,14]. In this paper, the functional form of the 
geodesies when A and r g are both non-zero are found, in terms of Riemann's 
hyperelliptic theta functions[Baker, 9,10]. 



3 Derivation of the Hyper-elliptic Integral 

Let a particle of mass m be constrained to move in the equitorial plane, 9 = 0, 
and let the action be given by [Landau, 15], 

S = -E t + L <j>+ S r (2) 

where E is the total energy of the particle and L is it's angular momentum. By 
the conservation of four momentum, 

dS dS 2 /q\ 
9 dx* d& ~ m ^ 
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we can insert S and g lJ from ([!]) and (Q)into (||) and obtain, 

_ r (E 2 - (1 - 1 A a r 2 - r g /r) (L 2 /r 2 + to 2 c 2 ))l 
" " L {l-\Kar 2 -r a /r) * (4) 

Differentiating SV w.r.t. L and using 



= const = (5) 



(6) 



ro r 2 (E 12 - (f - | A a r 2 - r g /r) (L 2 /r 2 + m 2 c 2 ))2 

This integral is a hyperelliptic integral of genus, g = 2. The problem is to 
invert this integral, i.e. express the upper limit r as a function of <fi >[J a ~ 
cobi, 16], [Gauss, 17], [Abel, 18], this is known as the inversion problem. 

4 The Solution to the Hyper-elliptic Integral of 
Arbitrary Genus 

The following will be a review of the account given by H.F. Baker[9,10] and 
the references contained in these books. The general solution to the inversion 
problem for hyperelliptic integrals is given by solving g of the 2 g + 1 equations, 
for the g variable places x\, . . . ,x g , 

V WU ' -A(b) (b- Xl )...(b-x g ) (7) 



The notation Baker uses for the generalized theta function is, 

$(u\u b - a ) = n m . m ,) = ~ m, i to') = e a 6(«; ~ to, ~ m') (8) 

where a is an arbitrary g x g symetrical matrix, since putting (^J) into (Q) gives, 
d 2 {u\u b ^) _ 9 2 («;i m,| to') 



7? 2 (u|it 6 '.°') e 2 (w; i fc, i k>) 
where to, to', /c and k' are integers. Riemann's theta functions are defined as, 

Q( v - 1 m _ TO M = J] e 2 k " m ')+ b («+3 ™') 2 +i T m m') Qq'\ 

V ' 2 ' 2 

with, S = X)ni°^-oo • ■ ■ Sn°^-oo- his&gXg matrix,in general non-symmetrical. 
b is a g x g symmetrical matrix. In (||), u denotes the g quantities, 

(11) 
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with, i = 1 . . . g. This is known as Abel's Theorem. Riemann's Normal Integral 
of the first kind is defined as, 

( x : l)i,g— 1 d&_ (12) 



(x, = X s - 1 + A l:9 _ 2 x 9 2 + ... + A i>0 (13) 



where 



and Ai g -i, . . . , Ai o are arbitrary constants. The denominator in (|I2| ) is given 

by, 

y 2 = 4 (x - ai) . . . (x - a g ) (x - a) . . . (x - c g ) (x - c) (14) 

where ax, . . . , a g , ci, . . . , c are branch places and x is a variable on the one com- 
plex dimensional Ricmann surface of genus, g, shown in figure 1. 
From figure(El) the period table, table 1, is derived. 



Table 1. Primativc Normal Pcriods[Brioschi,19] [Bolza,20] [Konigsberger,2I] 



Int 






a g 


Hi 




Pg 


x.a 

u{ 








w' n 




u 'ig 
















x, a 


UJgl 




Ugg 


u'gi 




<9 
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Now let, 



mi w r< i + ... + m g u r ^ g + m' rl ui' rl + ... + m' g u' rg 



(15) 



where m and m' are equal to +1 when integrating across a period loop from 
right to left, -1 when crossing from left to right, and when not crossing a 
period loop. Now define an integral of the first kind via, 



7T i vp a = h rA u{ a + h ra u x 2 ' a + ... + h r ^ g u x g ' a 
with r = l,2,...,g and where, 

2 h w = tt i 

b = 2hw' = -KiT 
so with these definitions table 1 becomes table 2, below 



(16) 

(17) 
(18) 



Table 2. Hyperclliptic Moduli 



Int 




a 2 






01 






x.a 


1 










n,i 






















x,a 










1 






T g,g 



with Riemann's condition for hyperelliptic integrals, [Brioschi, 13], [Bolza,14], 

r = r* (19) 

With this condition there are \ n (n + 1) + 1 different non-zero periods. With 
the above equations, Riemann's theta function becomes, 

0( w - 1 m , - m!) = £ e 2 v i v ( n+ i m ')+ l 71 T ( n +i ™') 2 +* w ™ (™+| m ') (20) 

where the quadratic forms are given by, 

v m ) = v n+i m' v = v\ m + . . . + v g n g + ^ vi m' 1 + . . . + ^ v g m' g (21) 



and similarly, 

T (n+- m 'f 



t n 2 + t n m! H — r m' 2 
4 

= (n,i n 2 + 2 ti j2 n\ n 2 + . . . + r SiS (n s ) 2 ) 

9 9 



+EE( 

s— 1 r — 1 



t t . s ) n r m s 



+ 4 (n,i Ki) 2 + 2 n, 2 mi m' 2 + . . . + r 9 , 9 K,) 2 )(22) 
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It has been proven by Riemann that for the integral of the firt kind, i.e. v^ ,a , 
when it has the period scheme as in table 2, then the imaginary part of, 

71,1 n\ + 2 Ti )2 fix n 2 + ■■■ + T g .g n 2 g (23) 

is positive for all integer values of ri\ and n 2 , with the exception of m = n 2 = 
0. Therefore the modulus of e 1 * T n is less than unity and the function (p0|) 
converges for all values of it's argument, v. The last thing to define in @) is the 
constant A(b), it is defined as, 



A{b) = (e ^- {x-a 1 )...{x-a g ) (x - a) . . .(x - c g )(x - c))\ x = b ) a (24) 



where, e is +1 when u b,a is an odd half period and -1 when it is an even half 
period. 

5 The Solution to the Hyperelliptic Integral of 
Genus Two 

The solution of the hyper-elliptic integral^, is given by(0),[Rosenhain,22],[Gopel,23],[Forsyth,24], 
[Baker, 9], with g = 2 we have 



<d 2 {vi,V2] \ mi, \ m 2 , \ m[, \ m' 2 ) 
Q 2 (vi,v 2 ; 5 h, i k 2 , 5 k[, \ k' 2 ) 

where by(|l6|)we have, [Baker, 9, 10], 



= A{b) (b-xi) (b-x 2 ) (25) 



x,a i x,a . 7 x,a 

IT I — /li,i U x +Al,2 u 2 

x i a j. x,a . j x,a 

■K I V 2 = «2,X % + «*2,2 U 2 (26) 

where /i^j are the elements of the following matrix, 

7T i / UJ 2 o — Wl,2 



h 

2 A V — W24 wi,i 
and, 

A = det(/i) ^ (27) 



Abel's Theorem(|Tl|) , becomes 



U x ' + % ' = Wi 



j, 2 -r u, 2 

these become with, 



x±.ai i 2:2^1 /no^ 
U + Wo = «2 (28) 



x 2 = a 2 (29) 



G 
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Figure 2: 
Genus Two Riemann Surface 



and with (12 



Ul 

rXl (A lfi + A hl x) dx 



y 

(A 2 , + A 2 ,i x) dx 

Ol V 

(30) 



where, 



j/ 2 = 4 + A 4 x 4 + A 3 x 3 + A 2 x 2 + Ai a: + A (31) 

These coefhents may be substituted into the general solution to the quintic 
equation [Drociuk,25] where the branch places, i.e. the roots of the quintic (|3l|), 
are determined in their functional form. Then (|3l|) maybe rewritten in the form, 

y 2 = 4 (x - ax) (x - a 2 ) (x - c\) (x - c 2 ) (x - c) (32) 

So the the Riemann Surface, figure 1, becomes a genus 2 Riemann surface figure 
2. 

From figurc(H) the period table, table 3, is derived. 



Table 3. Period Table for Genus Two Hyperelliptic Integrals 



Int 


OL\ 


a.2 


01 


02 


x.a 








U'12 


x.a 

u 2 


W21 


^22 


u'21 


w 22 
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With definitions (|2q), table three becomes table four, 
Table 4. Genus Two Hyperelliptic Moduli 



Int 




OL2 


01 


02 


it. a 
v l 


1 





n,i 


71,2 


x.a 
V 2 





1 


72,1 


7"2,2 



So the four non-zero periods are, 1 and 

1^2,2 ^11 — Wl,2 W 2 i 



Til 



72,2 



and by Riemann's condition(|l9|), 



7-1,2 = 



-U>2,1 ^2.2 • W 2 2 



-W 2 ,l tL>i 4 + Wl,l W 21 



k>2,2 &12 ~ w l,2 w 2,2 



= 7- 2: i 



(33) 



(34) 



The quadrupoly periodic Double Theta Function is given is from([20D , (|2 1|) , (|22D , 
B(wi, u 2 ; - mi, - m 2 , - m l7 - m 2 ) 



+mi)+u 2 (2 n 2 +m' 2 )) , 



ni— — oo ?i2 — — oo 



^(2 n 2 +"i2) 2 q'( 2 "1+mi) 2 ^(2 rii+roi) (2 n 2 +m' 2 ) x ^_ jjmj ni+ro 2 «2 x (j)" l i " l 'i+™2 ™ 2 



where, 



(35) 
(36) 



q = e* "-^ 
q' = ei l7T T2 ' 2 
r = e i 1 T i,2 



(37) 



which is equal to Forsyth's double theta function, [Forsyth, 24], with the two 
characteristics related by the following, 

mi = p 
m 2 = p' 



m, = a 



(38) 
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so that a r in Forsyth's paper is equal to 

a r = {i) a ' p ' +p " r^ 2 " +cr ') ( 2 m + si a ma ) (39) 

where n — n\ and m = n 2 are the summation indices. Forsyth gives a review of 
Rosenhain's theory, [Forsyth, 26] of the fifteen ratio's of the quadrupoly periodic 
theta function. From the list given in [Forsyth, 24], the first ratio is selected and 
is equal to, 



Q 2 (vi,v2\\ (1), I (1), I (1), I (0)) = / (oi - gg) , _ ) 

e a («i,«2;J (1),| (0),i (0),| (1)) V (a x - d) ( ai - c 2 ) (a x - c) [Xl ° } 

(40) 



where (|24|) (|25[) (p9[)and( 35 ) have been used. The double theta funcntions used in 
@ are 

oc oc 

g 13 = e »i (»i (2 ni + l)+^2 (2 n 2 )) ? (2 n 2 ) 2 

m=— 00 ri2 — — 00 

? /(2 ru + 1) 2 r (2 m + 1) (2 n 2 ) /-gm+ria { £ 4 jj 

where Forsyth's notation is adopted, 

9i3 = 0(vuv 2 ; |,|,|,0) (42) 
and in the denominator of (^Q), is 



00 00 

t i in 1 > I — ;.., r> n-,4- II) (1 

ni— — 00 n 2 — — 00 

/(2 m) 2 (2 ni) (2 n 2 + l) 



J 12 _ gTT i («i (2 ni)+u 2 (2 n 2 + l)) „l2 ! i 



^'(2 m) r (2 m) (2 n 2 + l) (43) 

0i2 = e(«i,« 2 ;i 0,0, i) (44) 

6 The Solution to the Hyper-elliptic Integral 

The solution to the hyper-elliptic integral(0)is, with xi = — , 



r = (r g B\ 2 y/ai - a 2 )/(c \[a\ - a 2 0\ 2 + V ( a i - d) ( a i - d) (ai - c) 6^ 3 ) 

(45) 

where #13 and #12 are given in ([4l|)and(fi"3|). The hyper-elliptic integral(|6])may 
be rewritten as, 

^ * (46) 
<„, 2/2 

and define a new variable, (/>' , given by 

r*i dx 6' 

7 - 2 ( 4? ) 
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so that the variables in the argument of the theta functions are, 



v x = ((/i M A lfi +h ia A 2 , ) (/>'- (h lt i A u +/i li2 A a ,i) 0)/(2 7T t) 

U2 = ((/l2,l Ai )0 + /l2,2 A 2 , ) </>' - (/l 2 ,l Ai,i + /l 2) 2 A 2) l) 0)/(2 7T i) (48) 



where h is defined in (pq)and(p7|). The constants are given by, 

Ai.o = = A 2 , (49) 

Ai,i=A 2l i = l (50) 

With this choice the equations of Baker and Forsyth are identical to these, in 
the case of g = 2. The roots of the quintic,ai, 02, c\, c 2 and c are given by setting 
the coeffients in ([H]) equal to, 

A = (4 A to 2 c 2 r*)/(3 L 2 ) (51) 

Ai = (52) 

A 2 = 4 (^ 2 c 2 L 2 + A/3 - m 2 c 2 ) r 2 (53) 

A 3 = (4 m 2 c 2 r 2 )/L 2 (54) 

A 4 = -4 (55) 

then substuting into the solution to the quintic [Drociuk, 25]. If one thinks 
in terms of the advancement of perhelia, i.e. Mercury, in the case A = 0, the 
orbits obey doubly periodic functions. In the case when A ^ they are equal 
to quadrupoly periodic functions of two variables <p and <f>' , equations ( |46| ) and 
(Et]). <t> is the physical angular coordinate defined in equation ([!]). <fi' exist 's 
because of Abel's theorem, equation (Bq). Simply taking 0' equal to zero cannot 
work for then equation (|27])will go to zero and the transformation (|2^) becomes 
singular. Therefore cj>' is an unknown parameter? 



7 The Reduction of the Accellerating Universe 
to the General Relativistic Universe 

When the Cosmic Force is set to zero and particles are subject to gravitation 
alone, two roots of the quintic, see [Drociuk, 25] or by inspection, go to zero. 
Then ti. 2 = t 2 .2 0, then q' — r = 1, and setting n — 0, gives 

§3 ^ _ fc sn 2 (vi) (56) 

^12 

Where sn is one of Jacobi's four elliptic functions [Forsyth, 26]. Choosing to let 
a-2 = c-2 => 0, then ([i"5|) becomes 

- k sn 2 (vi) = (at - x{)/y/ (ai - ci) (ai - c) (57) 
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W U ^1,1 




Figure 3: 
Genus One Riemann Surface 



where, 

a>i = rg/n 

ci = r g /r 3 

c = r g /r 5 

k = ^(ci-aiV^-ai) (58) 

which gives the Schwarzschild black hole solution without the cosmological con- 
stant, as a genus 1 elliptic function, which is the same as,[Bartlett,14], it is 

(---) = (---)sn"(v 1 ) (59) 
r n r 3 n 



Since two roots are zero, (|32|)becomes with 02 = C2 = 0, 

y 2 = 4 (x - a x ) [x - c x ) {x - c) (60) 

which gives rise to genus 9=1 functions whose Riemann Surface is, Figure 3 
and the period table reduces to, 



Table 5 Normal Primative Periods Genus One 



Int 






x.a 







then by (|26|), with = ^2,1 = ^2,2 = 0, gives 

/i = h\\ = 7r i/u>n 

Then table 4 reduces to, 



(61) 
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Tabic 6 Genus One Moduli 



Int 


ai 




x.a 


1 


Til 



Where, 



and 



ni = Ki)/(wu) = ( / dx/y)/( dx/y) 



-(— / A x ,i ete/y) 

^11 7 a 



which becomes with (46) 



« x = 0)/(a;ii 2) 



and choose, 



i4 M = w u (r 9 (l-l))Va 



(62) 
(63) 

(64) 
(65) 



8 Conclusion 

The geodesies of particles in the combined Schwarzschild black hole and a vac- 
uum energy density, are given by ratios of quadrupoly periodic theta functions 
of two variables. The four moduli of the theta functions are given by ratios 
of periods of the hyperelliptic integrals of the first kind. These these primitive 
normal periods may be obtained by integrating between the branch places on 
the Riemann surface of genus two. 

Acknowledgement: I would like to thank the University of British Columbia 
for the use of the computing facilities to produce this document. 
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